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In the paper [2j we studied the temporally inhomogeneous system of non-colliding Brow- 
nian motions and proved that multi-time correlation functions are generally given by the 
quaternion determinants in the sense of Dyson and Mehta. In this report we give another 
proof of the equivalent statement using Fredholm determinant and Fredholm pfaffian, and 
claim that the present system is a typical example of pfaffian processes. 
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1 Non-Colliding Brownian Motions 

By virtue of the Karlin-McGregor formula jSJ Ej, the transition density of the absorbing 
Brownian motion in a Weyl chamber 



jx = (xi, x 2 , ■ ■ ■ , x N ) G M N ; x\ < x 2 < • • ■ < XnJ, 



is given by 



/jv(t;x,y) = det 

l<t,j<N 



x,y6E<, t£[0,oo) 



-(x-y) 2 /2t 



'2nd 



where p t is the heat-kernel given by 

Pt(x,y) -- 

The integral 

Mv(t;x) = 

where dy = YliLi ^Vh gi ves the probability that a Brownian motion starting from x G 
does not hit the boundary of IR< up to time t > 0. 



/jv(t;x,y)dy, 



a JV 



For a given T > 0, we define 



g NT (s,yi;t,y) 



f N {t - s; x, y)A/jy(T - t; y) 
A^(T-s;x) 



for < s < t < T,x,y G R< . It can be regarded as the transition probability density from 
the state x G M< at time s to the state y G M< at time t. In jHJ it was shown that as 
|x| — > 0, g^(0,x;t,y) converges to 



AT 



;i-2) 



^, T (0, 0; t, y) = C(N, T, t)h N (y) JJp t (0, y»)A/jv(T - t, y) 
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where C(N,T,t) 



71 



N/2 



nf=i r(i/2) ) 7 , Ar(JV-i)/4^-AT(AT-i)/2^ and 



l<i<j<N 



det 



The N particle system of non-colliding Brownian motions X(t) all started from the origin at 
time 0, i.e. X(0) = = (0, 0, • ■ ■ , 0), and conditioned not to collide with each other in a time 
interval (0, T] is defined by the process associated with the transition probability density g^,T 
given by (jl.l)) and (jl.2j) . This process is temporally inhomogeneous and it was obtained as a 
diffusion scaling limit of the vicious walker model in Figure 1 illustrates the process X(t). 




Figure 1: Process X(t),i G [0,T], with X(0) = 0. 

Let 3C be the space of countable subset £ of R satisfying JJ(£ D -ft') < oo for any compact 
subset .fT. For x = (xi, X2, ■ ■ ■ , x n ) G l^J^ =1 IR £ , we denote {xj}" =1 G X simply by {x}. The 
diffusion process {X(t)}, t G [0, T] on the set X is defined by the transition probability density 
Q N>T (s,{^};t,{y}),0<s<t<T: 

r &v,T(s,x;*,y), ifa>0, tf{x} = |}{y} = 
fljv,T(s,{x};t,{y})= < 9N, T (0,0',t,y), if s = 0, {x} = {0}, |j{y} = N, 

y 0, otherwise, 

where x and y are the elements of ]R<. For the given time interval [0, T], we consider the 
M intermediate times to = < t\ < £2 < • • • < tM < tM+x — T. The multi-time transition 
density of the process {X(i)} is given by 

M 

QNtT (t , {x<°>}; t u {x (1) }; . . . ; t M+1 , {x( M+1 )}) = J] QN , T (t„ {>)}; t, +1 , {x^}). 

/i=0 
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Here we set t = 0, t M +i = T and {x (0) } = {0}. From (JHU) and (JOJ) we have 

(1-3) 0jv,t(O, {0}; t h {x«}; . . . ; t M+1 , {x^ 1 )}) 

= C(N,T, tl )h N (x«)sgn (h N (x^ 1 ))) 



JV M 



(a») _(a*+i) 



For xW G R£, 1 < n < M + 1, and iV' = 1, 2, . . . , N, we write x$ = (x[^,x^\ x ( $ 

For a sequence {-/V^KJj 1 of positive integers less than or equal to N, we define the multi-time 
correlation function by 

(1-4) Pn,t (fi, {x«}; f 2 , {xjg}; . . . ; t M+1 , 



n <^y, n 

S/V , T (0, {0}; t 1; {x«}; i 2 , {x£>}; . . . ; t M+1 , {xjf . 



Let Co(M) be the set of all continuous real functions with compact supports. For f = 
(A, h, ■ ■ ■ , /m+i) e C (K) A/+1 , and = (0i, 2 , • • • , M +i) G E M+1 , the multi-time character- 
istic function is defined for the process X(t) = (Xi(t), X 2 (t), • • • , Xjv(i)) as 



{M+l TV 
v^5>m£/m(*«,(*m)) 
jU=l i M =l 

where E N ^ T [ - ] denotes the expectation determined by Qn,t- Let 

Xli {x) = e^WrW - 1, 1 < n < M + 1, 
then by the definition of multi-time correlation function (jl.4|) . we have 

JV JV 



*jv,T(f; o) 



N 

E E- E 

JVi=0JV 2 =0 JV M+1 =0 

JVf+1 



dx 



dx 



(2) 
JV 2 



dx 



(M+l) 



■JV 



V M+1 



M+l 



(1.5) 



II II {x«}; f a , {xg>}; ■ ■ ■ ; t M+1 , {xg£?}) . 



2 Fredholm Pfaffian Representation of Characteristic 
Function 



Let 



z n,t[x\ 



1 

M 

A4 



Jtf+1 p M+l 

J R N(M+1) 



FT dxW det [M^Ocf K(0, )(1 + xi(a:f ))" 

/z=l 



l<i,3<N L J \ 



(M+l)x 
JV > 
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where Zjv,t[0] is obtained from Z n ^t[x] by setting Xn( x ) = f° r a U /•*• 



Here Mj(x) is an arbitrary polynomial of x with degree i in the form Mi(x) = biX % + • ■ • with 
a constant 6j for i = 0, 1,2, and thus h N {yi) = deti<jj<jv Mj_i(xj) /IlfeLi^c-i- Then 
()1.3|) gives 

( 2 1) iTr iv. a\ _ Z N,t[x] 

en 

By the well-known formula 



/ dx. det 


fa{Xj) 


det 


M X j) 


= det 


/ dx 4>i{x)(f)j{x) 


/ R w l<ij<JV 


l<i,j<N 




l<i,j<N 



for square integrable continuous functions fa, (pi, 1 < i < N, we have 



dy det 

l<i,j<N 



x 



/ TT ^ M.^^K^x^^i + xi^)) 

u=l ^ J 



where p T _ tM+1 {x,y) = p {x,y) = 6(x-y). 

For simplicity of expressions, we will assume that the number of particles N is even from 
now on. We use the formula 



dy det 

jv l<iJ<N 

< 



My j 



pf 



l<i,j<N 



dy / dysga(y - y)fa(y)fa(y) 



for integrable continuous functions fa, 1 < i < N. Here, for an integer n and an antisymmetric 
2n x 2n matrix A = ((%), the pfafnan is defined as 



Pf (A) = Pfl<,<i<2r 
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^ Sgn(a ) 0^(1)^(2) Q(t(3)(t(4) • • • «<T(2n-l)<r(2n), 



where the summation is extended over all permutations a of (1,2,..., 2n) with restriction 
<r(2A; — 1) < a(2k),k = 1,2, ... ,n. Since (Pf(A)) 2 = det A for any antisymmetric 2n x 2n 
matrix A, we have 



^tv,t[x] 



det 

l<i,j<N 



X 



*'/' /f '"" f ' " 

/ TT dxWiM i _ 1 (^)p tl (0, 2 ;W)(l + xi(x (1) )) 

M f . 
W+l ^ 



M f . 

x II + X,+i(^ +1) )) \pt- 

(1=1 { > 



*M+1 



(x^\y) 



det 

l<i,j<N 



(A ) tJ + (A 1 ) ij + (A 2 ) ij + (As) 



M+l 

E E 

^=1 l<//i<// 2 <---<W<M+l 



y dysgn(y-y) J dxM i _ 1 (x)p tl (Q,x)p T _ tl (x,y) 

y dy sgn(y-y) 

x / fjtfa^) / dx Mj_i(x)p tl (0, x)p tMi _ tl (x, x^) 
J& k=1 J 

x n (»^ ft) )\ +r tjx (ft i> +i) )}»(> ) )^ ( (^ w) .!/) 

k=l ^ ' 

x y dxMj_i(x)p tl (0,x)p T -ti(5,y), 



M+l 



Y dy d v s ^ n (y-y) 

l=l i< t i 1 < t i 2 <-<ne<M+l J J 
x / TT dx^ k) [ dxMj^WjHM^Ptn-tM^) 

Jul ^1 J 



e-i 

x 
k=i 

M+l M+l 



YY Y Y I dy I d v s ^ n (y-y) 

£=1 rn=l l<fj, 1 <fi 2 <--<(i e <M+ll<ui<U2<---<u m <M+l' 



X 



X 



f 

[ f[dx^ [ dxM i _ 1 {x)p tl {^x)p t ^ tl {x,x^) 
J& k=1 J 

k=i ^ > 

m „ 
n=l J 

- 1 r i 

X [I X n (^K +1 A(^^ K+l) ) \Xrn(x M )pT-uJx^\y), 



n=l 
m—1 



where we have used the Chapman-Kolmogorov equation for the heat-kernel 

Jdy Pt -,(x,y M ,z)= Pt ,-,( X ,z), 0<s<t<u, X ,yeR, 

and the fact that po(x, y) = \im t _> p t (x, y) = 5(x — y). 

We consider a vector space V with the ortho normal basis < x) \ which 

I ) fx=l,2,-,Af+l,xeR 

satisfy 

(ji,x\i/,y) = 5 fM/ 5(x - y), 

fi, v — 1, 2, • • • ,M + l,x, t/GK. We introduce the operators J, p, p+,p~ and x acting on V as 
follows 

(//,x|J|i/,y) = l( /1=v =Af+i)Sgn(y - x), 
(/i,x|p|z/,y) =p| t ,_ tM |(x,y), 
(2-2) = Pt„-t„(a;,j/)l(,*<„) = x), 

where is the indicator function: 1^) = 1 if is satisfied and = otherwise, and we 
will use the convention 

M+l . 

= / dyA(p,x;is,y)B(is,y;p,z) = {p,x\AB\p,z) 

for operators A, 5 with (/i, x|A|z/, y) = A(p, x; u, y), (p, x|-B|z/, y) = B(p, x; is, y). 

Consider another basis | \i) ; i = 1, 2, • • • j in V and we assume that transformation matrix 
between the two bases is given by 

(2-3) {i\l*,x) = (p,x\i) = J dyM i . 1 (y)p tl (0,y)p t ^t 1 (y, x), 

% = 1, 2, • • • , p = 1, 2, • • • , M + 1, a; G R. Then the quantity 
y dy M i _ 1 (y)p tl (0,y)p t ^ tl (y,x) 

ii— 1 p m „ 

/TEPm / 

m=i i<m</i 2 <-</im</i i=i 

m-l . . 

n XM fc (^ (Mfe) K fe+1 -^(^ (Mfe) ^ (Mfc+i) ) kj^Wwo^*), 



X 

fe=l 

for p — 1, 2, • • • , M + 1, x e K, i = 1, 2, • • •, can be written as 



(i\p,x) + ^(il^x^HfjL^x^lxp+l^x^) 



m>l 



■ ■ ■ (/i m _i, X^ m - l] \XP+ | A*m, ^ (Mm) ) (/Xm, X) 

(<||i,x> + X;<i|(XP+) m |^a;> 



m>l 



It is also expressed as (u,x\ — 

i-p-x 

It should be noted that the basis | \i) ; i — 1, 2, • • • | is in general not orthonormal. Here 
we introduce an operator 8 such that 

(2.4) (i\6\j) = {j\6\i) = 6 ij , z,j = 1,2,---. 
We will use the following convention 

oo 

m 3 )( 3 \B\^x) = Y,^ 3 Bf\x) 

for Aij = (i\A\j) = (i\fi,x)(fi,x\A\u,y)(u,y\j),B i j ' 1 \x) = (J\B\fi,x) = (J\u,y)(v,y\B\fi,x), 

but we will not write it as (i\AB\/j,, x), since ; i — 1, 2, • • • | is in general not a complete 

basis. By this basis, any operator A on V may have a semi-infinite matrix representation 
A = ( (il-Ah') ) . If the matrix A representing an operator A is invertible, we define the 

V ^ 1.2,.. 

operator A A such that its matrix representation is the inverse of A; 

(2.5) ((*\A A \j)).. =A-\ 
In other words, 

(i\A\j)(j\A A \k) = (i\A A \j)(j\A\k) = (i\8\k), i,k = 1,2,-... 
For any given operator A, the equality 

(i\A\j)(j\A A \k)(k\B\e) = (i\B\e). 

holds for arbitrary i, i — 1, 2, • • • and B. Then the equality 

(2.6) A\j)(j\A A \k)(k\ = l 

should be established for each A. We will use this equality later. 

Let Vn be a projection operator from the space spanned by ; i = 1,2, • • • j to its 

iV-dimensional subspace spanned by ; i — 1, 2, • • • , TV j, and thus 



(il^jvl^x) = (n,x\V N \i) = 



(i\fx,x), if 1 < i < iV, 
0, otherwise. 

Then we have the following expressions for (A*)ij, « = 0, 1, 2, 3, i, j = 1,2,---, N, 



(A )ij = (i\V N \n,x){ii,x\j\is,y){v,y\V N \j) 
= (i\V N JV N \j), 



7 



it 



'J 



(As) 



V 



1 



-x\fi,x)(fi,x\pJV N \j) 



N~ ~ 

V n x- — |//, x) (//, x\pJV N \j) 

i-p+x 

Vnx- —pJV N \j), 

i -p+x 

V N Jp\/j,,x)(/j,,x\x 
1 



V N Jpx 



i -p-x 



v N \j), 



V 



1 



-r^-xl/^, x\pJp\u, y) (v, y\x- —V N \ j) 

i-p-x 

^jvX- — 1^, s|pJp|z/, y) (v, y\x~ —V N \j) 

i - P+X 

Vnx-, — —pJpx- 



1 -p-x 



l-p+X i- P-X 



V N \j). 



That is, the matrices A a — ((A 

representations in the basis ||z) ; i = 1, 2, • • • | of the operators. 



, a = 0, 1, 2, 3, can be regarded as the matrix 

■ ' i,j=l,2,-,N 



Since (z N)T [0]\ = det 



l<ij'<JV 



(2.7) |*iv,T(f;0) 
By our notation (|2.5j) . 

since it satisfies the relation 



det 

l<i,j<N 



(A )ij\ , flUD gives 

5ij + {A^Ax)ij + (Ao 1 A 2 ) ij + (A^As)^ 



(A % = (t\V N (V N JV N ) A V N \j), 



(i\VN(V N JV N rV N \j)(j\V N JV N \k) 

N 



/](Aq 1 )ij(^o)jfc — l(l<t,fc<JV)^ifc 



*-{l<i,k<N) 

(?|7Wiv|j>. 



We will use the abbreviation 

A N = V N AV N 

for an operator A. Then it is easy to see that ()2.7)1 is written in the form 

2 



(2.8) 



det 

l<i,j<N 



<5y + (i|B|/i,x)(/i,x|C|j) 



det (i| 

l<i,j<N 



5 N + BC 
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where (i|B|/^, x) and (/i, x|C|j) are the two-dimensional row and column vectors, respectively, 
given by 



{i\V N J^V N \k){k\x T= ^-E\^x) -(i\V N J^V N \k)(k\Jpx\p,x 



A, 



(i|B|/i,x) 
<MC|j) 



-p+x 



\ 



-(i\V N J^V N \k) 



*{^Xi=^\v,y){v,y\pJpx\^,x) J 



{n,x\pJV N \j) 

i^x\ TZ ^r-*V N \j) 

The determinant ()2.8|) is written using the Fredholm determinant, 



Det(/i, x\ 



/ 2 + CB 



Det 



+(H,x\pJ\i) 

*(t\VNjNV N \3)(j\ XT =^\v,y) 
x(i\V N j£v N \j)(j\ XT 4^\v,y) 



~(p, x\pJ\i) (i\V N J^V N \j) (j\JpxW, y) 
-{^,x\0\i}{i\V N J^V N \j)(j\XTz^j\p, z) 

x( P ,z\pjpxW, y) 

(p,,x\v, y) 

+(^,x\ T 4—^\i)(i\V N J^V N \j)(j\JpxW, y) 



Det 



(fi,x\p,z)^ 
+(/j,,x\pj\i) 

x(i\V N J^V N \j){j\xj4^\p,z) 
x(i\V N J^V N \j)(j\ XT 4^\p,z) 



+ (p,^T=jZ±\i)(i\'PNJ N V N \j)(j\xT^-^\p,z) 
x(p,z\pjpx\v, y) 

-(P,x\pj\i) 

x(i\V N J^V N \j)(j\Jpx\p,z) 
+ (p,x\pJpx\p, z) 



{p,x\p, z) 

x(i\V N J^V N \j)(j\Jp X \p, z) 

(p,zW,y) -(p,z\pJpxW,y) 







(P,z\u, y) 



Det 



(fi,x\u,y) 
+ (p,x\pJ\i) 

x(i\r N j£v N \j}(j\ XT 4-zW,y) 



-(P,x\pj\i) 



x{i\V N j£v N \j)(j\Jpx\v,y) 

+(fi,x\pjpxW, y) 

(p,,x\v, y) 

x(i\V N J^V N \j)(j\JpxW,y) 



x(i\v N j^v N \j){j\xT^m\^y) 
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where I2 denotes the unit matrix with size 2. It is further rewritten as 
Det(fj,,x\ I J 2 + 



1 




i-p-x 



pj\i)(i\v N j£r N \j)(j\x ( -pJ\i)(i\v N J^v N \j)(j\jp+pjp)x 



\i)(i\V N j£v N \j)(j\ X 



\i)(i\V N j£v N \j)(j\jp X 



r 1 


0' 


) 





1 





W, y) 



Det( / u, x| 



l-p-x 



1 



+ 



pj\i)(i\v N j£v N \j)(j\x ( -pJKXipjvj^WOVp+^jx 



i)(i\V N J^V N \j)(j\ X 



i)(i\V N J^V N \j){j\Jp X 



Det(^,x| I 2 + 



pJ\i){i\V N J A V N \j){j\ -p+ -pJ\i)(i\V N J^V N \j){j\Jp + pJp' 



\t)(i\PNj£PN\j)(j\ 
Here we have used the facts that 

"1 

Det(/i, x\ 



i)(i\V N J^V N \j){j\Jp-p. 



and 



Det(/i, x| 



l-p-x 
1-p+X 



k y) 



W,y) = 1, 



1 

which are consequences of definitions ()2.2j) of the operators p + and ]3_ . Then we arrive at 



(2.9) 

where 
(2.10) 



jv,T(f ; <?) ^ = Det J 2 - y) + 



S»> v {x,y) P*' v {x,y) 
D^ u (x,y) S u ^(y,x) 



Xu{y) 



D^(x,y) = -(v,x\i)(i\V N {V N jV N ) A V N \j)(j\v,y), 
S^(x,y) = (^x\pJ\i){i[P N {V N JV N ) A V N \j){j\^y), 
I^(x,y) = -(^x\pJ\i)^\V N (V N JV N ) A V N \3)( 3 \jp\^y), 



and 
(2.11) 



S» v (x,y) = S" v (x,y)-(iJL,x\p + \v,y) 

(fi,x\pj\i)(i\V N {V N JV N ) A V N \])(]\v,y), if ii > v, 

-{fx,x\pJ\i){i\(^J A -V N {V N JV N ) A V N j\j){j\p,y), if //< 1/, 

10 



I^ u (x,y) + (n,x\pJp\v, y) 

(li,x\pj\i)(i\( J A -V N (V N JV N ) A V N )\j)(j\Jp\v,y), 



I^(x,y) 



where we have used the equality (|2.fi|) for A = J. 

In (121 Rains introduced Fredholm pfaffian, denoted here by PF, and proved a useful 
equality 

(2.12) |PF(J 2 + K)| = Det(J 2 + J^K), 
for any antisymmetric 2x2 matrix kernel K, where 

(2.13) J 2 = 



1 
-1 

Then ()2.9|) implies the Fredholm pfaffian representation of the multi-time characteristic func- 
tion „ 

D"' u (x,y) S u '"(y,x) 



(2.14) ^(f; 9) = PF J^Jix - y) + 



-S^ v (x,y) -iw{x,y) 



Xv{y) ■ 



3 Pfaffian Process 



Let 
(3.1) 



A^(x,y) 



D^ u (x,y) Sy(y,x) 
-S»> u (x,y) -I^ v (x,y) 



and construct 2 J^u^i x ^ SulV antisymmetric matrices 



a (.A) x (2) ... x (w+iA _ / A w( r to Jy) 



l<i<N^,l<j<N„,l<fi,v<M+l 



for N m = 1, 2, • • • , TV, 1 < m < M + 1. By the definition of Fredholm pfaffian jTH] and the 
equality (|2.14j) . we can establish the following statement. 

Theorem 1 The non-colliding system of Brownian motions X(£) is a pfaffian process in the 
sense that any multi-time correlation function is given by a pfaffian 



(3.2) 
(3.3) 



Pn,t (h, {x^}; t 2 , {x^}; . . . ; t M +i, { x NjJ^} 



Pf 



(M+l)^ 



Remark. The pfaffian processes considered here may be regarded as the continuous space- 
time version of the pfaffian point processes and pfaffian Schur processes studied by Sasamoto 
and Imamura ^Ej and by Borodin and Rains pQ . The processes studied in || are also pfaffian 
processes, since the 'quaternion determinantal expressions', in the sense of Dyson and Mehta 
[3 El Ell- of correlation functions are readily transformed to pfaffian expressions. 



11 



Let Hi(x) be the i-th Hermite polynomial 



HAx) = e* 2 ~ e~* 2 



d 

dx 



K (2xY-^ 



where [a] denotes the greatest integer not greater than a. The Hermite polynomials satisfy 
the orthogonal relations 

(3.4) / dxe- x2 Hi(x)Hj(x) = 2 i i\^8 ij , i, j = 0, 1, 2, • • • . 

Set 



t 1 (2T-t 1 ) I2T - t x 

ci = v^ — ' zi= v~^' 

and 

{2~ l c\Sij, if % is even, 

2- i <*{6 ij - 2{i - l)8i- 2j }, if i is odd. 

Now we specify polynomials {Mj(x)} as 
(3.6) M i {x)=b i z^Y^a ij H j (-\zl i = 0,1,2,- •• 

f 1-1/2 1 /^2\2i+l/2 

with h = |r [i/2] | , where r< = ~r(i + l/2)r(z + 1) ( ^ J .Set 

where 7^/2 denotes the unit matrix with size iV/2 and J 2 is given by (|2.13jh and let J be the 
semi-infinite matrix obtained as the N ^ oo limit of J^. By the orthogonality of Hermite 
polynomials ()3.4jl . we can show through ()2.3j) with the choice (|3.6|) that jHUZj 

(z|J|j) = (i\ii,x)(n,x\j\v,y)(v,y\j) = J ij: i,j = 1,2,---. 

Since Jjr = —1^ for any even N > 2, this implies that the matrix ( (i\ J\j) ) is invertible 

and 



(3.7) (i\J A \ 3 ) = =-J t 



(3.8) (i\V N (V N JV N ) A V N \j) = (i|7W A ?W) 

= -{J N )ij, if 1 < i,j < N, 



0, otherwise 



for i,j = 1,2, 
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If we write 

(fi,x\i) = biB^ x {x), 
(i\Jp\fi,x) = -(fi,x\pJ\i) = 

% = 1, 2, • • • , /i = 1, 2, • • • , M + 1, x G R, then the functions (|2.1()jl are written as 

JV/2-1 



(3.9) D^foy) = 

S^(x,y) = 

and Equations (j2.11|) become 



i=o 1 

N/2-1 

i=o n 

N/2-1 



^(s)^i(v)-^i(*W(v) 



>(")/ 



i=0 



(3.10) S^(x,y) 



N/2-1 



8=0 



^(*)^i(y)-^i(*)^(y) 



(*0 



?(")/ 



x 



i=N/2 



if yU > z/, 



if // < 



/"•"(x,y) 



x 

i=AT/2 * 



Theorem ^ with the expressions ()3.9j) . (|3.1()jl of the elements of matrix (j3.1j) is equivalent 
with Theorem 3 reported in although the latter was given in the form of quaternion 
determinant. The present argument will be generalized to discuss other non-colliding systems 
of diffusion particles reported in [TQ1 [TT] . 
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